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I. INTRODUCTION 

This is the sixth paper in a series of studies of Par- 
tial Quenching (PQ) in Chiral Perturbation Theory 
(%PT) Q, 0, U , where we now investigate the effects of 
including electromagnetic loop corrections in the theory. 
The motivation for using partial quenching in xPT (and 
indeed for studying xPT itself) comes from the fact that 
even though quantum chromodynamics (QCD) over time 
has become the generally accepted theory of the strong 
interaction, it has still proven difficult to use this theory 
to derive low-energy hadronic observables such as masses 
and decay constants. 

An alternative approach is to use Lattice QCD simula- 
tions to this end. Computational limitations have how- 
ever hindered such simulations for light particles since 
they can propagate over large distances, requiring very 
large lattice sizes. Because of this, most simulations have 
so far been performed with heavier quark masses than 
those of the physical world. 

On the other hand, Chiral Perturbation Theory pro- 
vides a theoretically correct description of the low-energy 
properties of QCD, and can be used to extrapolate the 
results of lattice simulations down to the masses of the 
physical regime of QCD. In particular, one can use lat- 
tice simulations to determine the low-energy constants of 
%PT by fitting x?T calculations to corresponding lattice 
simulations, thereby getting estimates of hadronic low- 
energy observables. 

One problem with this approach is that reliable predic- 
tions from xPT require that one keeps the quark masses 
fairly small, and so far it has proven difficult to reach 
the chiral regime in the lattice simulations. However, 
progress is being made on this front. A complementary 
approach is therefore to use partial quenching where one 
introduces a separate quark mass for the calculations of 
closed quark loops, referred to as sea quarks, compared to 
the quark lines which are connected to external sources, 
referred to as valence quarks. This has the advantage 
over full QCD calculations that results with more val- 
ues of the valence quark masses can be obtained with a 



smaller number of values of sea quark masses, which is 
useful since varying the latter is computationally more 
expensive. 

Unquenched QCD may be recovered from partially 
quenched QCD (PQQCD) by taking the limit of equal 
sea and valence quark masses, and therefore it follows 
that QCD and PQQCD are continuously connected by 
the variation of sea-quark masses. This means that, in 
contrast to a fully quenched theory where the effects of 
the closed quark loops are neglected altogether, one can 
relate partially quenched QCD simulations to the un- 
quenched physical observables of the real world. 

Chiral Perturbation Theory has also been extended to 
include both quenching and partial quenching 0,|HE1;0- 
The formulation of Partially Quenched x?T (PQxPT) 
is such that the dependence on the quark masses is ex- 
plicit, and thus the limit of equal sea and valence quark 
masses can also be considered for PQxPT. This allows 
for determination of the physically relevant LECs of xPT 
by fits of partially quenched xPT (PQxPT) to partially 
quenched lattice simulations (PQQCD), see e.g. the dis- 
cussion in In particular, the LECs of xPT, which are 
of physical significance, can be obtained directly from 
those of PQxPT. More detailed discussions of this and 
references to earlier work can be found in the papers of 
Sharpe and Shoresh |8|,|9j- The calculations in this paper 
have been performed in three-flavor PQxPT without the 
$0 HI degree of freedom. In our earlier work with Timo 
Lahde Cfl, IH El 13 H we have calculated masses and 
decay constants for the charged, or off-diagonal, mesons 
to next-to-next-to-leading-order in PQxPT. However, in 
order to compare with the experimental values at high 
precision one needs to take electromagnetic effects into 
account as well. 

Electromagnetic corrections in xPT have a long his- 
tory. The lowest order (LO) was done by Dashen 0. 
The first corrections to this were worked out in Ref. |l6j| . 
That larg e corrections might appear was pointed out in 
Refs. flTl ITsL Il9j . where these corrections appear both 
from chiral logarithms and effects persisting at large 
N c . The work in pure xPT was started by Urech|20j 
and by Neufeld and Rupertsberger [2l]]. The NLO ex- 
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pressions for the masses were calculated in both papers 
and the decay constants in the second. More recent 
work on estimating the relevant LECs can be found in 
Refs. IH il HI |H IH Note that there are sev- 
eral subtleties involved in defining electromagnetic cor- 
rections at low energies as discussed in Refs. jU H3, [H] 
with increasing levels of detail. There are also first ex- 
ploratory lattice QCD calculations |29l l.'S(l l3lL l32| . 

In this paper we present the extension of PQxPT to 
include dynamical photons to NLO. In addition we calcu- 
late the NLO expressions for the masses of the charged, 
or offdiagonal, mesons with virtual photon-loop correc- 
tions for all possible degrees of degeneracy in the valence- 
and sea-quark masses. We do the same for the decay 
constants, i.e. we determine the 0(e 2 ) and 0{e 2 p 2 ) cor- 
rections in PQxPT to all these quantities. 

We point out that the two phenomenologically rele- 
vant quantities AM 2 and AF for differences of masses 
and decay constants, can be determined from partially 
quenched lattice calculations where the photons are only 
coupled to the valence quarks. This has the important 
consequence that these quantities can be calculated in 
lattice QCD with gluon configurations generated with- 
out dynamical photons. 

The paper is organised as follows. First we present 
the technical background and notation already present 
in the earlier work [H GH H 13 El in Sect. |TT| We 
also present the results for the needed loop integrals 
there. The extension to dynamical photons is discussed 
in Sect. IIIII Here we give the Lagrangians needed, as well 
as the subtractions needed to obtain finite results. The 
analytical expressions for the masses and decay constants 
are given in Sect. HVI and discussed in Sect. [V] Some il- 
lustrative numerical results are give in Sect. IVII and we 
recapitulate the main conclusions in Sect. IV ill 



II. PQxPT, TECHNICAL OVERVIEW 

Here we give a short overview of the technical aspects 
of PQxPT. A more thorough discussion of the technical 
aspects of PQxPT calculations to NLO (without pho- 
tons) can be found in Refs. Our earlier papers, in 
particular > a l so contain overviews of the NLO techni- 
cal and notational details, but the focus there is mainly 
on the NNLO aspects relevant for those papers. Lectures 
on standard x^T can be found in |33| . 

The mechanism which gives different masses to sea 
quarks and valence quarks in PQxPT is introduced by 
adding explicit sea quarks, as well as unphysical bosonic 
ghost quarks. The bosonic quarks are needed to cancel all 
effects of closed loop contributions from valence quarks. 
This cancellation happens if the masses of the bosonic 
quarks are set equal to the masses of the valence quarks. 
The symmetry group of PQxPT is essentially given by 
the graded group 

G = SU(n val + n scei \n va \) L x 5f7(n va ] + n soa |n va i)_R . (1) 



where n va \ denotes the number of valence quarks and n sea 
the number of sea quarks in the theory. The number of 
bosonic quarks is by necessity equal to the number of 
valence quarks. The PQ analog to the field matrix U in 
ordinary xPT is given by 



U = exp 



(iy/2$/F 



(2) 



The matrix $ is now a graded matrix, which in terms 
of a sub-matrix notation for the flavor structure can be 
written as 



$ = 



qvqv 



qsqv 



qBqv 



qvqs 



qsqs 



qsqs 



qvqs 



qsqs 



qsqB 



\ 



(3) 



The brackets denote matrices of the form 

u a Ub u a dp u a Sb 
't.,% = I d a u b d a d_ b d a s b 
s a Ub s a d b s a s b 



(4) 



where we have used three quark flavors u, d and s and the 
labels V, S and B in the sub-matrices stand for valence, 
sea and bosonic quarks, respectively. In general, the size 
of each sub-matrix depends on the exact number of quark 
flavors used, but for this paper all blocks in Eq. © are 
3x3 blocks 

The quarks qv, qs and their respective antiquarks are 
fcrmions, while the quarks qs and their antiquarks are 
bosons. Each sub-matrix in Eq. therefore consists of 
either fermionic or bosonic fields only. This construction 
means that $ satisfies the usual rules for cyclicity under 
trace and determinant products, provided that we per- 
form the corresponding sypersymmetric operations in- 
stead. For the trace, we must instead take the supertrace, 
defined by 



Str 



A B 
C D 



TvA-TiD. 



(5) 



where A, D, denote block matrices with commuting ele- 
ments and B, C denote block matrices with anticommut- 
ing (fermionic) elements. 

This also has the very useful consequence that the La- 
grangian structure of PQxPT is the same as for rt-flavor 
xPT, provided that the traces of matrix products in those 
Lagrangians are replaced by supertraces. A detailed dis- 
cussion about the Lagrangians and LECs for the different 
versions of xPT and PQxPT without the <E>o can be found 
in [13|. This correspondence between n-flavor xPT and 
PQxPT also holds for the divergence structure when re- 
placing n with the number of sea-quarks. The same also 
holds for the extension to electromagnetism in the next 
section. 

In the following, the different quark masses are iden- 
tified by the flavor indices i = 1, . . . , 9, rather than by 
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the indices it, d, s and V,S,B of Eqs. I@J and 0. The 
results are expressed in terms of the quark masses m q via 
the quantities \i = 2£?o m qi , where B is related to the 
quark-anti-quark vacuum expectation value in the chi- 
ral limit. Thus XiiX2,X3: belong to the valence sector, 
X4jX5jX6 to the sea sector, and X7iX8,X9 to the ghost 
sector. Since we set the quark masses of the ghost sec- 
tor equal to the quark masses of the valence sector, the 
masses X7?X8 7X9 do not appear in the final analytical 
results. 



A. Loop Integrals and Notation 

The expressions for the NLO masses and decay con- 
stants of the charged pseudoscalar mesons depend on 
several loop integrals. The renormalized contributions 
from these integrals are written in terms of the functions 



Mx) = 

5(x,x;0) = 

B'{Xi,X-,x) = 

#i(x 7 ,x;x) = 

B[(x7,x;x) = 



-7Ti 6 xlog(x/M 2 ), 
-Trie (1 + log(x/M 2 )) , 

7T 16 (1 - l0g(x/M 2 )) +0(X 7 ), 

•0(*r). 



^ f 1 + \ log(x/x 7 ; 



7T16 



log(x//i 2 ) + G(x 7 ), 



tt 16 /(2x) + G(X 7 ), 



(6) 



where fj, denotes the renormalization scale and Trie = 
l/(167r 2 ). The argument x 7 is a small photon mass intro- 
duced to regulate infrared divergences. The prime indi- 
cates a derivative with respect to the momentum squared 
in the loop integral. 

The quantities d va i and rf S ea are used to indicate the 
number of nondegenerate quark masses in the valence 
sector and the sea sector respectively. For d va i = 1, one 
has xi = X2 = X4, while d va \ = 2 means xi = X2 ^ X4- 
d V ai = 3 is not needed for this paper. Similarly, d sca _ — 1 
means X4 = X5 = Xe, d sea = 2 means X4 = X5 ^ Xe and 
for d sea = 3 all the sea quark masses are nondegenerate, 
such that X4 7^ X5 7^ Xe- 

The lowest order neutral pion and eta meson masses 
in the sea quark sector show up at several places in the 
analytical results. They are denoted by x-n and Xn and 
are given by the relations 

2 / 

Xtt + Xv = 3 U4 + X5 + Xe) , 

XttXt? = 3 (X4X5 + X5X6 + X4Xe) , (7) 

which have no polynomial solution for d sc& = 3, but for 
d S oa = 2 one has Xtt = X4 and Xn = (X4 + 2xe)/ 3 - For 
d sca = 1 this simplifies further into x-n — Xv = Xi- The 
neutral meson propagators in PQxPT generate certain 
reoccuring combinations of the sea and valence quark 
masses. An overview of this can be found in ,13]. The 
relevant quark-mass combinations for this paper can be 



expressed in terms of the general quantities R* b defined 
by 



abc 



R 



abed 



R. 



abedefg 



Xa ~ Xb, 
Xa - Xb 
Xa Xc 

(Xa - Xb)(Xa ~ Xc) 

Xa ~ Xd 
(Xa - Xb)(Xa - Xc)(Xa ~ Xd) 
(Xa ~ Xe)(Xa ~ Xf)iXa ~ Xg) ' 



(8) 



and so on. For the case of d sca = 3, the needed combina- 
tions are 



R 



jkl 

T?d 



R 



i4567T?7 ) 



n ijkl 



Rjkk + RjU 



t" ^6x77 ^-tttiti ^-kttii 



(9) 



For the case of c? soa = 2, corresponding combinations are 



R 



jk 

Rt 



H6jk > 

z 

i46r/ s 



and for cL 



1, one has 



(10) 



R. 



pa pz 



1. 



(11) 



For certain sums and differences of quark-masses, or 
electric quark charges, we introduce shorthand notation 
given by 



Xi 



q, :! 



~ i = 4,5,6 

Si - Qj, 

i = 4,5,6 



(12) 



The quark charges are expressed in terms of the unit 
charge e. qij is the charge of a meson with flavor quantum 
numbers of quarks qiqj. 

The summation conventions from Ref. |l3j | have as well 
been implemented where possible. In short, they are as 
follows: 

• If the index s is present, the entire term is to be 
summed over all sea quark indices. 

• If the index q is present in a term, there will always 
be an index p and the resulting sum is over the pairs 
of valence indices. If only p is present, the sum is 
just over the valence indices. If we choose valence 
quarks of type 1 and 3, this becomes summing over 
(Pi q) = (1; 3) and (p, q) — (3, 1) or if only p is 
present, the sum is over p = 1 and p = 3. 
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• If the index m is present, there will always be an 
index n and the corresponding sum is over the pairs 
(m,n) = (77,77) and (m,n) = {t),tt). If only the 
index n is present, then the term is to be summed 
over the x?r and Xv masses. 

III. VIRTUAL PHOTONS 

In xPT, photons are included as external vector fields 
through a charge matrix Q for the three light quarks 
and through the covariant derivative D M j2£j- Introduc- 
ing photons in PQxPT is completely analogous, provided 
that traces are replaced by supertraces, in the following 
denoted by (• • • ), and that we use the n-flavor expressions 
for the Lagrangians. The covariant derivative includes 
the photon field through 

D^U = d„U -ir^U + iUl^ (13) 

with 

r M = v f _ l + eQ R A^ + 

In = Vn + eQ L A^ - a M . (14) 

For the meson masses, we set w M = = and for the 
decay constants = 0. The charge matrix Ql,r is the 
natural generalisation of the SU (3) charge matrix in [20| . 
e is the absolute value of the electron charge for physical 
quantities but is a free parameter in the lattice calcula- 
tions. The L/R notation refers to the symmetry prop- 
erties assigned to the Q's during the construction of the 
allowed Lagrangian terms |20j| , but in the usual physical 
picture the charge matrix is a constant matrix, and one 
has Ql = Qr = Q-, where Q is a diagonal matrix given 
I'.v 

Q = diag(<ji, . . . ,q g ). (15) 

However, for the notation used below, the distinction be- 
tween the two types is still needed. Furthermore, one 
would normally set q\ =2/3 and 92 = 93 = — 1/3 to agree 
with ordinary photon-included xPT for the real world, 
but for greater generality we have kept the g,'s free in 
the analytical results in this paper. In ordinary xPT one 
requires the charge matrix to be traceless. For PQxPT, 
Q is a graded matrix where we set qi = q%, <fe = q%, and 
99 = 93- This, together with the earlier requirement on 
the masses insures that the closed valence quark loops 
with photons coupled to them also cancel against the 
corresponding ghost quark loops. Therefore the PQxPT 
requirement (Q) = becomes 

?4 + 95 + <?6 = 0. (16) 

Finally, the quark masses are present through the matrix 



For convenience, the Lagrangians below will be written 
in terms of the field matrix 

u ee exp (i$/(V2Fo)) , (18) 

which is related to U through u — \/TJ . We also introduce 
the quantities 

Uft = i{u^(d li — ir lt )u — u(d li — il li )w}, 
X± = w f X ut ±MX ( u, 

Q L = uQlm 

Qr = vJQru 

V m Q l = uDpQiu 1 

V^Qr = u+D^Qru, (19) 

where Ft and Fr denote the field strengths of the exter- 
nal fields I and r, such that 

F ,xy = guy -d u l>* -i[l>*,l v ], 

= d"r u -d u r ti -i[r fi ,r v ], (20) 

and the covariant derivatives of Ql, Qr are defined by 

D»Q L = d"Qt-i[l^,Qt], 

D»Qr = O^Qr — i [r**, Qr] ■ (21) 

The quantities in Eq. (I19|) have a well-defined and simpler 
(as elaborated in |13ll34l |) behaviour under the symmetry 
transformations needed for the construction of the La- 
grangians. In this notation, the lowest order Lagrangian 
has the form 

£2 = - ~ V" - ~A(d^) 2 
F 2 

+ e 2 C(Q L Q R ), (22) 

where F^ is the field strength tensor of the photon field 
A^, with F^ u — d^A v — d^A^. Furthermore, A is the 
gauge fixing parameter, here set to A = 1, and e is the 
electric unit charge. We will also use the notation Ze = 
C/Fq. The lowest order Lagrangian contains terms of 
0(p 2 ) and 0(e 2 ). 

For £4, the result is as well analogous to [2(j, except 
that the terms presented there are for SU(3). For the 
n-flavor case needed in PQxPT, one has two additional 
LECs due to the fact that the Cayley-Hamilton relations 
needed for the derivation of £ 4 ^ only are true for the 
SU(3) case. We split the NLO Lagrangian into the purely 
strong part of 0(p 4 ) and the part including electromag- 
netic interactions up to 0(e 2 p 2 ), and thus write 



X = diag(xi, • ■ • ,Xs), Xi = 2B om qi . (17) 



£4 — £g4 + 



(23) 
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The strong part is given by 



valence and sea-quark masses equal is 



Csi = LiXi + contact terms 



i=0 

= L Q {U^U V U^U„) + L 1 (U^U^) 2 + L 2 (M M Ii")(tl p H,;> 

+ L 3 ({u^u^f) + L 4 (u'V) (x+) + L h (u"u M x+) 
+ U ( X+ ) 2 + L 7 < X _) 2 + ^ ( X 2 + + X 2 -) 

P^ . „. \ i LlO i t 2 f 2\ 



#i (F 2 +F 2 )+i7 2 ( xx t), 



(24) 



where the Li and £/i are the partially quenched LECs for 
the case with three sea quark flavors. 

The electromagnetic part to 0(e 2 p 2 ) is 



C S 2E2 = e 2 F 2 | > , KfQ\ + K&QU + K&Q' W 




with 



Qi 

Qi 
Qi 
Qi 
Qt 
QI 

Qt 

QI 

Qg 

Qio 
Q'u 

Ql2 

Ql3 
Ql4 
Ql8 
Ql9 



- (Ql + Qi) 

(QlQr) (u M u") 
-(Ql?v) (Qlu" 

(Ql<v) (fiRU") 

((2^ + 2^/) 
((QlQr + QrQl)<V^> 

\{Ql + Ql) (x+> 

(QlQr) (x+> 
«Q L + Q R )x+) 
((QlQr + QrQl)x+) 
((QrQl - QlQr)*-) 

i([v M Q R ,Q R ] u"- [v m Q l ,Ql 

(V^QlV^Ql + V^QrV^Qr) 

+ Qrw m Qrw m ) 
(QlUuQru m ) . 



It") 



(25) 



For the Kf and Qf, we follow the numbering convention 
introduced by Urech HJ, but Kf 5 ,K^, Kf 7 are of 0(e 4 ) 
and are not needed here. The new terms, needed for the 
partially quenched case are thus named K E S and K E 9 . 

The relation with constants Ki of Urech when setting 



I<2 

K 5 
K 6 
K, 



= K? + K 



18 ) 



= Ki 



Ik 



1!) ' 



K E - K E 



= Ki 



K 



19 i 

K E — 0K E 
A 5 z/v 18 ' 

K E - K E 
^6 ^19' 



= Kf 



i = 7, 



,14. 



(26) 



The extra subtractions needed can be derived from the 
divergences of the n-fiavour case. We write 

1 



Kf 



(e c M ) 



-2c 



K 



Er 



167r 2 e 



(27) 



(28) 



with the dimension of space-time d — 4 — 2e and 

C =-i(ln(47r)+r'(l) + l) . 

The equivalent subtractions needed for the Li can be 
found in Ref. 0, 35| . We have derived the values of the ki 
from the ro-flavor results given in the appendix of Ref. j3|J 
after correcting an obvious misprint and rewriting the 
terms in our minimal basis. The ki are give in Tab.[fl 



i 




i 


ki 


1 





9 


1 

8 


2 
3 
4 




-z 


10 
11 
12 


1 3 ry 

~8 _ 4 Z 
1 

16 
1 

8 


5 


3 
8 


13 





6 




14 





7 





18 


3 
8 


8 
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TABLE I: The values of the subtraction constants ki of 
PQ X PT. 



A. Propagators and LO masses 

The propagators for the supersymmetric formulation 
of PQxPT can be found in Ref. [9|. For calculational 
reasons, they have here been translated from the Eu- 
clidean formalism into Minkowski space. The charged 
propagators are given by [9j 



iGUk) 



k 2 -M$ Aj +i 



(i + 3) 



(29) 



where Mq ^ denotes the lowest order mass of the meson 
tJ for i ^ j and the signature vector e,- due to the graded 
structure is defined as 



+1 for j = 1, . . . , 6 
-1 for j = 7,8,9. 



(30) 
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For PQxPT without electromagnetic interactions, the 
LO mass M^- is simply Xij = (Xi + Xj)/%- Electro- 
magnetic interactions modify the lowest order mass, and 
the new LO mass, which in the analytical results below 
is denoted as Xe,ij> can be read off from the 0(<I> 2 ) terms 
of the lowest order Lagrangian. It is given by 



Xe 



Xi 



2Ce 2 
1 o 



(<li -QjT 



(31) 



This is the mass that appears in the charged propagator. 

The neutral propagators can have a double-pole struc- 
ture and are more complicated. However, since they are 
charge-neutral, the lowest order propagator is unaffected 
by the inclusion of electromagnetic corrections in the the- 
ory. Explicit expressions for the lowest order neutral 
propagators can be found in p| . See also @ . 



IV. ANALYTICAL RESULTS AT NLO 



physical masses, are defined by the position of the pole 
in Eq. (J33J), 



■^phys 



E(M 2 



physi 



Xi), 



(34) 



where the expression for the self-energy £ is written as 
a string of terms denoting the 1PI diagrams of progres- 
sively higher order. The contributions start at NLO, and 
thus 



£(M 2 hys , Xi ) = S 4 (M 2 , Xi) + 0(p e , eV), (35) 

Note that we have used the lowest order mass Mq instead 
of Mpjj ys in £4 since the diagrams in that term are already 
of 0(p 4 ,e 2 p 2 ). The Feynman diagrams that contribute 
to T*4(MQ,Xi) with electromagnetic corrections included 
are shown in Fig. ^ 



The analytical expressions for the masses and decay 
constants are fairly short, and very similar in form. 
Therefore it suffices to give the results for the cases with 
rfsoa = 3 only. They can also be downloaded from [37I . 
Expressions for the cases with d sca = 2 and e? sca = 1 can 
easily be derived by taking the appropriate limits, i.e. 
X5 -> Xi for 4ca = 2 and X5,X6 -> Xi f° r 4ea = 1. It 
should be noted, however, that for the degenerate cases, 
all sums are still over the full set of indices, and further- 
more, since we only take limits of the masses, the charges 
qi are never affected by such limits. 



A. Masses 

The corrections to the lowest order mass of a charged 
pseudoscalar meson is obtained by calculating the self- 
energy corrections to the propagator in the interacting 
theory, usually written in terms of the Fourier transform 
of the two-point function 

iA(p)= [ d i xe i r"{n\T[9(x) ji $(0)i j ]\n), (32) 



where $y = qiqj denotes any of the off-diagonal mesons 
in the valence sector of PQxPT, and f2 denotes the vac- 
uum of the interac ting theory. The propagator resums as 
a geometric series |38j . giving 



iA(p) 



Ml 



S(P 2 ,X.) 



(33) 



where M§ denotes the lowest order mass of the meson 
which is being considered, and x% m ^ denotes the de- 
pendence of the self-energy on all the lowest order meson 
masses. The quantity £(p 2 ,Xi) receives contributions 
from the one-particle-irreducible (1PI) diagrams. The 





FIG. 1: The Feynman diagrams that contribute to T,4(Mq , Xi) 
with electromagnetic corrections included. A filled circle de- 
notes a vertex from the £2 Lagrangian and an open square 
denotes a vertex from the £4 Lagrangian. A straight line is a 
pseudoscalar meson and a wiggly line is a photon. 



We present the physical mass on the form 



Xe,: 



j(4)vi 
r 



+ 0(p 6 ,eV), 



(36) 



where Xe,ij — Mq is the lowest order mass, defined in 
Eq. J3TJ. The superscripts (v) and (s) indicate the val- 
ues of d va \ and d sea , respectively. It should also be noted 
that the results are given in terms of the lowest order de- 
cay constant F and the lowest order masses, since these 
are the fundamental inputs in PQxPT. To the accuracy 
we are working with here they can be replaced by the 
physical masses in the NLO correction. 

The NLO contribution to the charged pseudoscalar me- 
son mass with electromagnetic corrections is for d va i = 1 



and gL 



3 found to be 



<5( 4 > 13 = [A8L r 6 -2ALl] X iXi + mi-8Ll]xl 

-48e 2 F 2 Z E L r 4 q 2 12 xi ~ 16e 2 F 2 Z E L r 5 q 2 12X i_ 
-e 2 F 2 [l2K^ r +\2K^ r -l2Kf r -l2Ki r ] Q 2 Xi 
-e 2 F 2 [<<M<-4<] q 2 pXl 
+12e 2 F 2 K£ r q 2 2Xl 
+8e 2 F 2 [K^+Kf{] q 2 2Xl 
-e 2 F 2 [8^44iff;] 9192X1 

-l/3A( X m)KnXi ~ l/3A( X i)Rtxi 
+e 2 F 2 A( X i)q 2 12 
+2e 2 F 2 Z E A( X is)q 2 12 
-l/3B( X i,xi,0)R d lX i 
+4e 2 F 2 B(x 7 ,Xi,Xi)'7?2Xi 
-4e 2 i?' 2 Bi(x 7 ,Xi > Xi)8i2Xi. (37) 
For d va i = 2 and <i sca = 3 one has 

5^ 23 = [48L r 6 -24V 4 \xiXi3 + [l6L r s -8L r 5 }x 2 13 

~48e 2 Z E F 2 L r 4 q 2 3X i - We 2 Z E F 2 L r 5 q 2 3X i3 
-e 2 F 2 [UK^+UKt-UK^-UKt] Q2X13 



2 F 2 [AK, 



Er 



Erl 



<? P X13 



+e 2 F 2 [4K 9 Er + 4Jf5r] g 2 x P 
+ 12e 2 F 2 ^g 2 3 xi 
+8e 2 F 2 [^fo'+Xfr] «? 2 3 Xi3 
-e 2 F 2 [8Kf 8 r + 4Kf 9 r ] qi q 3X i 3 
-l/3A( Xm )K 13 Xi3 ~ l/3A( Xp )RP q7rvX 
+e 2 F 2 A( Xl3 )q 2 13 
+2e 2 Z E F 2 A(xis)qi s qi3 
-2e 2 Z E F 2 A( X 3s)q3sqi3 
+4e 2 F 2 5(x 7 , Xi3, Xi 3 )9 2 3 Xi3 
-4e 2 F 2 J Bi(x 7 ,Xi3,Xi3)'7 2 3Xi3- 

B. Decay Constants 
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(38) 



The decay constants F a of the pseudoscalar mesons are 
defined through 



(0\AZ(0)\Mp))=iV2p' l F a , 



(39) 



in terms of the axial current operator A£(0). In the fol- 
lowing the flavor index a has been suppressed for simplic- 
ity. The Feynman diagrams that contribute to the axial 
current operator at NLO, are shown in Fig. [21 

Diagrams of 0(p 4 ,e 2 p 2 ) also contribute to Eq. Ij39(l 
through the wave function renormalization factor \/~Z, 
since the expression for the decay constant of a meson is 



-ftphys 



F« 



F A {M 2 



phys 



,X,) + 0(p b ,eV) 



(40) 





FIG. 2: The Feynman diagrams that contribute to the axial 
current operator at NLO. A filled circle denotes a vertex from 
the £.2 Lagrangian and an open square denotes a vertex from 
the d Lagrangian. A straight line is a pseudoscalar meson, a 
wiggly line is a photon and a dashed line is an axial current. 



In Eq. 140|) , the subscripts of the matrix elements F in- 
dicate the chiral order and the lowest order contribution 
F 2 has been identified with Fq. The wave function renor- 
malization is given in terms of the self-energy diagrams 
by 



1 



9S(P 2 ,X,) 



dp 2 



(41) 



which becomes, when expanded such that all contribu- 
tions up to 0(p 4 , e 2 p 2 ) are taken into account, 



£' 

3 = 1 + T 



E' 



dH A (p 2 ,Xi 



dp 2 



(42) 
(43) 



The quantity £4 denotes the one-particle-irreducible di- 
agrams to 0(p 4 , e 2 p 2 ). Combining these expressions, the 
decay constant at NLO is then 



F, 



phys 



F + F 4 ( Xi ) 
3£ 4 (p 2 ,X*) 



F 



2 dp 2 



0{p 



6 e 2 p 4 



(44) 



A At 



Again it is sufficient to use the lowest order mass Mq 
in £4 since the diagrams in that term are already of 
0(p 4 , e 2 p 2 ). The analytical results for the decay constant 
are below given in the form 



Fphys — Fq 



J"(4)vs 

l+ T -E?r+0{ P \eV) 



(45) 



As for the meson masses, the superscripts (v) and (s) 
indicate the values of d V ai and d soa , respectively. 

The NLO contribution to the decay constant for a 
charged pseudoscalar meson with electromagnetic correc- 



tions is for d va i = 1 and d s , 



3 found to be 
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dependent combinations appear: 



(4)13 _ 



K, 



2^f 9 r ] Qi q 2 



+2e 2 F 2 [Kf + 
+2 e 2 F 2 Kf/g 2 2 
+e 2 F 2 [4Kg + 

+2e 2 F*B'(x 7 , Xl ,xi)q 2 12 Xi 
-e 2 F$B 1 (x~ ( ,Xi,Xi)q 2 i2 
~2e 2 F^B[ (x 7 , Xi > Xi )q 2 2 Xi ■ 



(46) 



For d 



val 



2 and cL, 



3 the result is 



f(4)23 



= +12^X1 + 4^X13 
+6e 2 F 2 [K^ + Kf]Q 2 
+2e 2 F 2 [Kf + *f "] 
+2e 2 F 2 ^ (Z 2 3 
+e 2 ^ 2 [4^f 8 '- + 2Kf/] 9l g 3 

-l/12l(x m )^„ 13 
+^(X P ) [l/6i?^ - l/12i$ 

+ l/4A(Xe,p S ) 

-1/125(x p ,Xp,0)^ 
+2e 2 F 2 B'(x 7 :Xi3,Xi3)<Z 2 3Xi3 
-e 2 F 2 i?i(x 7 ,Xi3,Xi3)'? 2 3 
-2e 2 F 2 i?i(x 7 ,Xi3,Xi3)'??3Xi3- 



(47) 



The term containing -R"„ nl3 is somewhat tricky to take 
the limit to the simpler mass cases. The form needed for 
the simpler cases can be found in Ref. or m H3- 



V. DISCUSSION OF THE ANALYTICAL 
RESULTS 

Our analytical results are finite. The renormalization 
obtained from the n-flavour divergences using the ar- 
guments presented above and in our earlier work, did 
cancel those from the loop diagrams. In addition, they 
agree with earlier PQxPT results when the electromag- 
netic parts are removed as well as with the known results 
for electromagnetic corrections when removing the par- 
tial quenching. We have used the definition of the decay 
constant with the axial current. This definition has an 
infrared divergence as can be seen also in our result. We 
have regulated that divergence with a photon mass x 7 - 
This is the definition which was used in Ref. |2j as well. 
How to relate this to measurable quantities can be found 
in Ref. H|. 

Which combinations of the new LECs can now be de- 
termined from lattice calculations? In the masses 5 in- 



Y 2 
Y 3 
Yi 
Y 



= K 



= K, 



Er 
1 

Er 



9 



K 



Er 
10 • 



-Kt-Kt 



2K 

K, 



Er 
5 
Er 



2Ki 



2K 

-2K 



Er 
10 
Er 
18 



2Kf{ 

K Er 
19 ' 



(48) 



These can be determined by varying the charges and 
quark masses separately. It should be noted that the 
sea quark charges only have a dependence via Y\Q 2 X\3 
with undetermined LECs. 

The decay constants depend on the combinations 



Y 6 = K 
Y 7 = K, 



Er 
1 

5 + A 6 



Er 



K 



Er 
12 ' 



(49) 



as well as on I4. It should be noted that the sea quark 
charges only have a dependence via Y$Q 2 with undeter- 
mined LECs. 

There is in addition dependence on the sea-quark 
charges in the chiral logarithms, but this dependence is 
predicted at NLO. 

The individual masses and decay constants depend on 
the sea quark charges. But since the overall dependence 
on the sea quark charges appearing with unknown LECs 
is simple we can easily make combinations where this 
disappears. We use here the notation 



M 2 (xi,x3,gi,?3) 



(50) 



to denote the mass of the meson with valence masses xi 
and X2 and valence charges q\ and g 3 . The quantity 

AM 2 = Af 2 (xi,X3,gi,93)-M 2 (xi,X3,93,<Z3) 

-M 2 (xi,xi,qi,q3) + a/ 2 (xi,xi,<73,93X 51 ) 

is especially useful. Only the electromagnetic corrections 
survive and the only dependence on the sea-quark charges 
is in some of the chiral logarithms. Since these contri- 
butions are independent of the LECs, they can be sub- 
tracted before making fits with lattice simulations, and 
hence do not present any problem in this respect. The 
quantity in Eq. (|51(l is also directly relevant for the vio- 
lation of Dashen's theorem 0, Il8j |. 



AM 2 D = 



m 2 K0 ) 



(ml+ -m 2 ). 



(52) 



Dashen's theorem states that the electromagnetic part of 
AMp vanishes. AM 2 becomes the electromagnetic part 
of AM 2 ) up to some very small electromagnetic corrrec- 
tions to the 7r° mass in the isospin limit. 

Similarly, differences of decay constants of particles 
containg valence quarks with the same charges have 
no dependence on the sea quark charges with unknown 
LECs. In particular this true for the difference of the 
pion and kaon decay constant. Wc define 



F(xi,X3,qi,q3,{qsea}) 



(53) 
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to be the decay constant of a meson with valence masses 
and charges as for M 2 above and sea-quark charges 
{q sea }- The quantity 



AF = F( X i,X3,gi,93,{<Z S ea})- J F 1 (xi,X3,0,0,{0}) 

-F(Xl, Xl> Ql, 93, {Qsea}) + F(xi, xu o, o, {0}) 
/*b , (54) 



is an example of this. It gives the relative electromag- 
netic corrections to the difference of kaon and pion decay 
constants. In fact, AF is independent of all the Kf r . 



VI. NUMERICAL RESULTS 

The whole purpose of this work is that our formulas 
can be used by the lattice QCD community to perform 
their fits. We therefore only present some representative 
numerical results. For the L\ we use the values deter- 
mined in the NNLO order fit of Ref. HjJ, called fit 10. 
For the extra electromagnetic parameters we use the es- 
timates of Ref. 22] . There are four combinations of the 
Kf r estimated there. We simply choose a series of Kf r 
values that reproduces the combinations estimated there 
and set all others to zero. The nonzero values we have 
chosen for illustration are 



K? r = 2.85-10" 



K 



Er 
10 

C 



4.0 ■ 1Q- 

ZeFq = 



-Er 
'9 
-Er 



K£ r = 1.3 • 10" 



iff/ = -1.25 • 10" 
4.2 • 10~ 5 GeV 5 , 



(55) 



at a subtraction scale \x = 770 MeV. Earlier estimates of 
C are in Refs. jHEI. 

The numerics we present here uses q± = q^ — 2/3 and 
<h = 93 = 95 = q& = —1/3 and a value of e determined 
from the measured fine structure constant a. We also 
only quote the electromagnetic part by subtracting the 
same result with e = 0. 

The lowest order correction to the meson masses van- 
ishes for those with zero total charge. For charged mesons 
it is equal to 



M 2 LO 



1.00 • 10~ 3 GeV 2 



(56) 



This should be compared to the physical mass difference 



1.3 • 10" 3 GeV 2 



(57) 



There is no electromagnetic correction to the decay con- 
stants at lowest order. 

Below we use for convenience the terminology tt for a 
meson with both valence masses equal to xi an d K for 
a meson with valence masses equal to xi an d X3 respec- 
tively. The charge label is + for valence quark charges 
2/3 and —1/3 and for valence quark charges —1/3 and 
— 1/3, i.e. electrically neutral. 



We first quote the electromagnetic corrections for xi 
Xa = X5 and X3 = Xe witn \/X~i = 135 MeV and y/xls 
495 MeV. 



^K+NLO 
M l°NLO 



= 0.45 • 10~ 3 GeV 2 

= 1.52 ■ 10~ 3 GeV 2 

= -2 • 10~ 7 GeV 2 , 

= -3 ■ 10~ 6 GeV 2 . 



This leads to a value of 



AM 2 = 1.07- 10" 3 GeV 2 



(58) 



(59) 



In agreement with the large violation of Dashen's theo- 
rem seen in Ref. |2^| since we used their estimate for the 
constants and similar values for the other inputs. The 
electromagnetic corrections for the decay constants with 
a photon mass x-y — (10 MeV) 2 are 



leading to 



F*+nlo/F = 0.0039. 
Fk+nlo/Fo = 0.0056, 



AF = 0.0017. 



(60) 



(61) 



The above results are for the unqucnched case. To 
show the effects of partial quenching we plot the quan- 
tities AM 2 and AF with input values as above and 
Xi = Xb an d X3 = Xe = 0.5 GeV 2 as a function of Xi an d 
Xi in Figs. 01 and ^ 



AM 2 , Fit 10, p 4 +e 2 p 2 




0.05 0.1 0.1 5 0.2 0.25 0.3 0.35 0.4 0.45 0.5 

/2i 



Xi 



FIG. 3: The quantity AM 2 of Ea.t5Tl the difference of electro- 
magnetic contributions to meson masses between kaons and 
pions as a function of the input lowest order masses xi an d 
Xi- The scale is logarithmic and contour lines are drawn at 
AM 2 = 0.00005, 0.0001, 0.0002, 0.0005, 0.001. 
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VII. CONCLUSIONS 



AF, Fit 10, p 4 +e 2 p 2 




0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.4 0.45 0.5 
Xi [GeV 2 ] 



In this paper we have shown how to include electro- 
magnetic corrections in partially quenched Chiral Per- 
turbation Theory. We have then used this formalism to 
compute the electromagnetic corrections to masses and 
decay constants of the charged or off-diagonal mesons 
to NLO in PQxPT. We also presented some illustrative 
numerical results. 



We have shown that for several phenomenologically in- 
teresting quantities the relevant LECs can be computed 
using quenched photons, i.e. they can be computed with 
the photons only coupling to the valence quarks. 



FIG. 4: The quantity AF of Eq. El the relative dif- 
ference of electromagnetic contributions to meson decay 
constants between kaons and pions as a function of 
the input lowest order masses xi an< l Xi- The scale 
is logarithmic and contour lines are drawn at AF = 
0.00005, 0.0001, 0.0002, 0.0005, 0.001, 0.002. 



The dependence on the sea quark mass is rather small 
in these differences. It cancels to a large extent. In Fig. El 
we show the electromagnetic contribution to the squared 
mass of the pion and kaon as a function of X4 — X5 with 
Xi = 0.1 GeV 2 and the other inputs as above. 



0.002 



0.0015 - 



> 

8 0.001 



0.0005 



It 

K + 



Acknowledgments 



1 1 1 1 1 1 1 1 1 1 1 

0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.4 0.45 0.5 

X 4 [GeV 2 ] 

FIG. 5: The NLO electromagnetic correction as a function of 
the sea quark mass X4 f° r the case of tv + and K + . See text 
for the other inputs. 
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